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Abstract. We report on Bloch-Zener oscillations of an ultracold Fermi gas in a tunable honey- 
comb lattice. The quasi-momentum distribution of the atoms is measured after sequentially 
passing through two Dirac points. We observe a double-peak feature in the transferred frac- 
tion to the second band, both as a function of the band gap at the Dirac points and the 
quasi-momentum of the trajectory. Our results are in good agreement with a simple ana- 
lytical model based on two successive Landau- Zener transitions. Owing to the variation of 
the potential gradient over the cloud size, coherent Stiickelberg oscillations are not visible 
in our measurements. This effect of the harmonic confinement is confirmed by a numerical 
simulation of the dynamics of a trapped 2D system. 

1 Introduction 

The transport properties of a quantum particle in a periodic potential depend crucially on the band 
structure of the system. A prime example is graphene, where remarkable transport behaviour results 
from the presence of two Dirac points in the band structure of its honeycomb lattice [1]. At a Dirac 
point, two energy bands intersect linearly and a quantum particle behaves as a relativistic Dirac fermion. 
Dirac points appear naturally in a number of condensed matter systems besides graphene, including nodal 
points in d-wdive superconductors or the surface states of three-dimensional topological insulators [2] . The 
transport properties associated to them can also be studied in artificially engineered systems using, for 
example, electrons in molecular graphene [3], microwave fields in meta-material structures [4] or ultracold 
atoms in optical lattices. 

Quantum gases in optical lattices have recently emerged as an attractive system for exploring Dirac 
points in a highly tunable environment. By loading a Bose-Einstein condensate into the excited bands 
of a one dimensional bichromatic optical lattice, the one-dimensional analogue of a Dirac point was 
realized and relativistic effects such as Zitterbewegung [5] and Klein tunnelling [6] were observed. Using a 
degenerate Fermi gas trapped in a two-dimensional honeycomb optical lattice, tunable Dirac cones were 
realized and the merging of Dirac points could be explored for the first time [7] . 

Another key feature of ultracold atoms is the possibility to perform transport studies in a dissipation- 
free environment, and to hence explore regimes which are difficult to access in solid state samples. For 
example, in the presence of a constant force, Bloch oscillations -the oscillatory motion of a particle 
in a periodic structure- can be observed [8]. In the presence of Dirac points, Bloch-Zener oscillations, 
a characteristic sequence of Bloch oscillations and Landau- Zener transitions between the two crossing 
energy bands, appear instead [9]. Their study allows for a precise characterization of the underlying band 
structure, giving access to properties such as the energy splitting between the bands [5] and the position 
of the Dirac points [7]. 
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In this paper, we present a systematic study of Bloch-Zener oscillations of an ultracold Fermi gas 
loaded into a tunable honeycomb lattice, focusing on quasi- momentum trajectories for which the two Dirac 
points are successively passed during the Bloch cycle. In the presence of a small energy gap at these band 
crossings each of them acts as an atomic beam splitter, where partial Landau-Zener tunnelling between 
the two energy bands occurs. Therefore, in this configuration the atomic wave-packet is successively 
splitted and then recombined during the Bloch cycle. Within the atomic cloud and when scanning system 
parameters we observe different behaviors depending on whether the probability of transfer at each band 
crossing close to a Dirac point is below or above 1/2. 

After presenting our experimental setup in Section _2, we explore the effect of the double Landau-Zener 
transition in Section 3. The experimental sequence is described in detail in Section 3.1. The fraction of 
atoms transferred to the second energy band is recorded as a function of the transverse quasi-momentum 
and the lattice geometry in Sections 3.2 and 3.3, and as a function of the sublattice offset in Section 3.4. 
We find good agreement with a simple analytical model based on two successive Landau-Zener tunnelling 
events. Despite the fact that the splitting process at the Dirac points is expected to be coherent, no 
interference fringes are observed, an aspect which is studied in detail in Section 3.5. We conclude in 
Section 4 



2 Experimental setup 

Recently, optical lattices with complex geometries have been implemented in various experimental se- 
tups [10-15]. The tunable optical lattice used in this setup is created in the same manner as described 
in Ref. [7J, combining in the x-y plane three retro-reflected laser beams of wavelength A = 1064 nm (see 
Fig. la). The two beams X and Y are arranged at an angle of 90.0(1)° with respect to each other and 
have the same frequency and polarization. They hence create a chequerboard lattice with a distance of 
A/ a/2 between potential minima. An additional beam X, propagating in the same direction as X but 
with a frequency detuning 5 ^ 0.4 GHz, gives rise to an independent standing wave of spacing A/2. 




Fig. 1. Experimental scheme, a Three retro-reflected laser beams give rise to an optical potential of honeycomb 
geometry. The lattice has a unit cell consisting of two sites A and B. b The corresponding 1^* and 2^^ Brillouin 
zones are shown, with the two Dirac points depicted as orange circles. Here qb = 27r/A denotes the Bloch wave 
vector, c Band structure of a honeycomb lattice for typical parameters. 



The overall potential is given by 
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V{x, y) = -V-j^cos^{kx + 6/2) - Vx cos^{kx) 

—Vy cos^ {ky) — 2a\/VxVY cos{kx) cos{ky) cos (1) 

where V^, Vx and Vy denote the lattice depths corresponding to each single laser beam (see Appendix B), 
k = 27r/A and the measured visibility of the interference pattern a is 0.90(5). The phase cp is stabilized 
to 0.00(3)7r, whilst can be continuously adjusted by tuning 5 and is set to values close to tt. The beams 
additionally cause a weak harmonic confinement for the atoms in all spatial directions, as described in 
Appendix B. 

The geometry of the lattice can be changed by adjusting the intensities of the three laser beams. In 
the following we focus on the honeycomb configuration, shown in Fig. la. The lattice has a unit cell 
consisting of two sites, labelled A and B. It gives rise to a band structure where the two lowest bands are 
connected by two Dirac points, but are well separated from higher energy bands. Since the angle between 
the primitive lattice vectors is set to 90°, the first Brillouin zone (B.Z.) has a square shape, and the Dirac 
points are located inside it (see Figs, lb and c). In the Qx direction, the refiection symmetry of the lattice 
fixes the Dirac points to the Qx — line. Their position in the Qy direction as well as the linear slope 
of the band structure in the vicinity of the Dirac points is determined by the strength of the horizontal 
and vertical tunnel couplings. The tunnelling can be controlled through the intensities of the laser beams 
forming the lattice. For example, increasing will mainly reduce the tunnelling in the x direction. 

Additionally, tuning the phase away from tt creates an energy offset Aab between sites A and B 
and hence breaks the inversion symmetry of the lattice. This causes a coupling between the previously 
orthogonal levels crossing at the Dirac points, and thus opens up a gap A there, which is proportional 
to Aab- 



3 Double transfer through Dirac points 
3.1 Experimental sequence 

We explore the band structure of our system using Bloch-Zener oscillations. The starting point of each 
measurement, as described in Ref. [7], is an ultracold cloud of N 50,000 fermionic ^^K atoms in 
the rriF = —9/2 Zeeman level of the Fhf = 9/2 hyperfine manifold. The gas is loaded into the lowest 
band of the tunable lattice. We set Vx/Er = 0.28(1) and Vy/Er = 1.8(1) unless otherwise stated. 
Here Er = /i^/2mA^, m denotes the mass of ^^K and h is the Planck constant. Along the third spatial 
direction z the atoms are only confined by a weak harmonic trap. 

We apply a force with characteristic energy Eb/Ji = FX/2h = 95(1) Hz to the atoms along the 
y direction using a magnetic field gradient, leading to Bloch oscillations [8] with a period of Tb = h/ Eb- 
After 11.5 ms (corresponding to roughly one Bloch oscillation cycle) the force is turned off. The optical 
lattice is then linearly ramped down in 500 ps and the gas is allowed to expand freely for 15 ms time of 
fiight before an absorption image is taken. This band-mapping technique [16-18] is used to measure the z- 
integrated quasi-momentum distribution of the atoms in the different energy bands. The final absorption 
image hence gives information about where in the Brillouin zone interband transfers have occurred. This 
can be seen in Fig. 2a, where for t ^ Tb atoms appear in the 2"^ B.Z. and are missing at the corresponding 
quasi-momenta in the 1^* B.Z. 



3.2 Landau-Zener transitions 

During one Bloch cycle the cloud probes the dispersion relation of the two lowest bands along the vertical 
Qy direction, which contains two avoided linear band crossings, see Fig. 2b. For Qx = 0^ the band crossings 
coincide with the two Dirac points. Neglecting the harmonic trap, one can treat trajectories for different Qx 
as independent. We attribute a transfer probability ^5 to a single band crossing along a (^^^ -trajectory, see 
Fig. 2_b. Neglecting the phase evolution and resulting interference of the two bands between the crossings 
(see Section 3.5), the probability for ending up in the higher band after passing the second crossing can 
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a t = o t==T^ 



b 




Fig. 2. Probing the Dirac points, a Absorption images of the quasi-momentum distribution of the atomic 
cloud after preparation (t = 0) and after performing one Bloch oscillation cycle (t ^ Tb). Colour scale: white 
indicates a high density of atoms, black indicates no atoms, b Cut through the band structure at qx = 0. During 
one Bloch cycle, the atoms pass through two band crossings. For qx = the band crossings coincide with the two 
Dirac points. Neglecting the phase evolution (see Section 3.2), each band crossing has a transfer probability ^s- 



be obtained by simply multiplying the probabilities for passage at the two crossings and summing up the 
two possible paths: 

= (1 - Cs) + (1 - Cs) = 2^5 (1 - Cs) • (2) 

The total transfer therefore has a maximum value of 1/2 for a Landau-Zener transition probability of 
= 1/2. 

Depending on the lattice parameters we either observe a single or a double slit of missing atoms in 
the lowest Brillouin zone after one full Bloch cycle, see Fig. 3a. The appearance of a double-peak feature 
in the transfer fraction £,d{Qx) is a direct consequence of the double transition through the two band 
crossings: for the lattice parameters on the left panel of Fig. 3a, ^5 < 1/2 at the central Qx = line, 
whereas ^5 > 1/2 for the right panel due to the different band structure, as will be detailed below. The 
single transfer probability ^5 generally decreases for quasi-momenta away from the center due to the 
increased gap at the crossing. Therefore, the maximum total transfer (i.e. = 1/2 <^ ^5 = 1/2) is 
reached at a finite value of Qx for the latter case. This gives rise to the double-peak structure. It appears 
due to the transfer probability being close to 1 at the two linear band crossings at the Dirac points. 

For a quantitative treatment, we perform a line sum of the optical density (OD) in the 1®* B.Z. (blue 
area in the figure). The initial cloud profile obtained from a fit to the quasi-momentum distribution after 
a Bloch cycle for a set of lattice parameters not giving rise to Dirac points is shown in purple. From these 
two profiles, the experimental transfer fraction (,d{Qx) (red curve below) can be obtained. As can be seen 
in the figure, the transfer is peaked ai Qx = for the left situation, whilst we observe a double-peak (with 
a dip at Qx = 0) for the situation on the right. 

The exact value of ^5 depends on the gap at the band crossing and the slope of the dispersion relation 
in the y direction far away from the crossing. To evaluate the transfer probability ^5 we use a simplified 
effective Hamiltonian well suited for describing the dispersion relation close to the two Dirac points, an 
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Vx,x,Y = [2.7(2),0.13(1),1.8(1)] Er 
9=1.01 3(1 )n 



Vx,x 



= [4.3(3),0.13(1),1.8(1)] Er 
9=1 .01 3(1 )n 
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Fig. 3. Transfer for different quasi-momentum-trajectories. a Transfer for two sets of lattice parameters. 
Left: Maximum transfer is observed for a central trajectory. Right: Maximum transfer happens to the left and 
right of the central trajectory. The plots depict the line sums along qy of the measured optical densities as a 
function of Qx (blue area) and the expected line sums without transfer (purple area). From this the transfer 
^d(Qcc) is calculated and compared to the prediction of the analytical model (red and black curve), b Numerical 
simulation of a trapped 2D system with N = 256 atoms: quasi-momentum distribution during one Bloch cycle 
for two exemplary situations. The population of the 2^^ B.Z. is plotted for two different -trajectories (position 
indicated by the coloured arrows). For the simulation parameters see Appendix A. 



approach previously worked out in detail in Refs. [19,20], 



H 



2m* 



A, 



2m* 



(3) 



Here = 



2/7,2 



denotes the merging gap, m* the effective mass along qy at the saddle point between 



the two Dirac points and and Cy the slopes of the dispersion relation along the different directions at 
the Dirac point, with h = h/27r being the reduced Planck constant. Their values are directly related to 
the lattice parameters [20]. Denoting the gap at the Dirac points with Z\, the transfer probability is then 



exp 



.^x^x 



ZiV4 



(4) 



The theoretical prediction of this simple analytical model is shown in the black curve in the bottom panel 
of Fig. 3a. The main features of the experimental momentum-resolved transfer fractions obtained from 
the atomic distribution in the lowest band are captured by theory. Deviations are possibly due to the 
finite resolution of the band-mapping technique as well as the uncertainties in the calibration of V^. 

The simple picture of two independent Landau-Zener transitions is confirmed when performing a 
numerical time-evolution of a trapped 2D lattice system in a tight binding limit, see Fig. 3b (details of 
the simulation see Appendix A). The transfer fraction at = increases with time on the left panel as 
the atoms pass the two Dirac points. On the right, however, after having passed through the first crossing, 
the transferred fraction decreases again. In contrast, the transferred fraction increases monotonously in 
both cases for quasi- moment a away from the center. 
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qx qx Dist. to B.Z. corner [qe] 

Fig. 4. Transfer for different lattice parameters, a Experimental (left) and theoretical (right) quasi- 
momentum distribution ^d{Qx) in the 1^* B.Z. after one full Bloch cycle integrated along the qy direction for 
different values of V^. The experimental data is the average of three consecutive measurements. The theoreti- 
cal quasi-momentum distributions are calculated under the assumption of a fixed ^^-distribution for the initial 
atomic cloud, which is obtained by a fit to the density profile at Vx/Er = 2.5 in the left panel. At this point 
the band structure does not contain Dirac points, b Extracted position of maximum transfer along qx using the 
distribution in the 2^^ band. Values and error bars denote the mean and standard deviation of three consecutive 
measurements The solid line shows the theoretical expectation of the maximum transfer position using the simple 
analytical model without free parameters. The experimental data was taken with set to 1.000(l)7r. 

3.3 Changing the slope of the Dirac cones 

As clearly visible in Fig. 3, the transfer depends on the precise lattice parameters. Since these values 
are freely tunable in our experiment, the dependency on the lattice structure can be systematically 
investigated. Fig. 4a shows such a scan versus along with the theoretical expectation using the 
simple analytical model. For V-^/Er < 3.4 there are no Dirac points in the band structure and thus no 
significant transfer is observed (no missing atoms). For V^/En > 3.4 two Dirac points are present in the 
band structure, leading to the formation of a progressively more pronounced double-peak feature as Vj^ 
gets larger. This is caused by the overall increased single transfer probability ^5 for larger due to the 
deformation of the Dirac cones, shifting the points of maximum transfer further apart. The position 
along Qx in quasi-momentum space for maximum transfer can in fact be obtained by extracting the peak 
position of the atomic distribution in the 2"^ B.Z. after taking a line sum along Qy. The results of this 
procedure are shown in Fig. 4b, which are in good agreement with the calculated position obtained from 
the simple analytical model. 



3.4 Opening a gap at the Dirac points 

The presence of a transition regime from a single- to a double-peak feature in the momentum resolved 
transfer S,d{Qx) depends on the gap at the Dirac points, which is controlled by the energy offset between 
neighbouring sites. For a vanishing gap only double-peak features appear, as the single transfer in the 
center is always close to 1. This is the case for the data presented in Fig. 4, as 6> ~ tt for this dataset. In 
contrast, an overly large gap merely leads to single-peak profiles with a very low overall transfer fraction, 
as the transfer in the center is already far below 1/2. Only if the gap is set to an intermediate value, 
a transition from a single to a double-peak profile occurs for increasing V^, which is the case for the 
parameters of the measurements in Fig. 3. 

From this simple picture we expect the total 3D averaged transfer fraction ^ to exhibit a dip as a 
function of the gap A at the Dirac point, since a large fraction of the atoms is located around = 0. A 
scan of the transfer fraction versus A for three different lattice configurations is shown in Fig. 5, where the 
transfer ^ for the entire atomic cloud was obtained by comparing the number of atoms after one full Bloch 
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Fig. 5. Scanning the gap at the Dirac point. The total transfer fraction ^ is measured for three different 
lattice depths versus the phase 0, which directly controls the energy offset between neighbouring sites. The 
gap at the Dirac point is obtained from a band structure calculation [7]. Values and error bars denote the mean 
and standard deviation of five consecutive measurements. The solid lines show the theoretical predictions from 
the simple analytical model without fitting parameters including the integration over Qx. 

cycle in the two lowest bands [7]. We observe a clear double-peak structure, which reduces for smaller 
values of V^. This is caused by the decreased slope Cy and increased slope Cx of the dispersion relation 
close to where the two Dirac points merge and annihilate, leading to an overall reduction of ^5 and thus 
to a less pronounced double-peak feature. The main features of the experimental results agree with the 
predictions of the simple analytical model using an integration over the entire cloud (see Appendix C). 
To calibrate the laser detuning S for which = tt we used Bloch oscillations along the x direction [7] . The 
slight shift of the symmetry axis for the data presented here is smaller than an estimate for a possible 
systematic error^. 

The comparison to theory allows to give an upper estimate for a possible small residual gap at the 
Dirac points. For the dataset where Vj^ = 3.6 Er a gap of 140 Hz would already lead to a vanishing of 
the double-peak structure. From this we conclude that, if there is any residual small gap at the Dirac 
points, for example caused by the finite size of the system, it is significantly smaller than 140 Hz. This is 
about a factor of 30 smaller than the bandwidth. 



3.5 Sttickelberg interference 

A full description of the system studied here would necessarily take coherence into account. The two 
sequential band crossings along one g^a^ -trajectory would therefore effectively realize a Stiickelberg inter- 
ferometer [21,22] with an associated dynamical phase 0. This phase depends both on the energy difference 
between the upper and the lower path (see Fig. 2b) and the time spent in between the two crossings, 
and therefore on the lattice geometry and the applied force. In the simple analytical model used so far 
this phase has been neglected. In fact, it can be taken into account by multiplying the transfer fraction 
with a correction factor 2cos^(^/2 + ^d), where the dynamical phase (j) and transfer phase (j)d are defined 
in [20] . Fig. 6 shows a comparison of the total 3D averaged transfer ^ versus for the simple incoherent 
model (blueTine) and the extended coherent model including this correction factor (purple dotted line). 
The experimental data agrees very well with the incoherent model, whereas the oscillatory behaviour as 
predicted by the extended model cannot be observed. 

For a correct treatment of the Stiickelberg interference, the variation of the effective force F due 
to the harmonic confinement over the cloud size has to be taken into account. Using an estimate for 

^ The estimate is based on an comparison to an independent calibration of the — n point using Raman-Nath 
diffraction on a ^^Rb Bose-Einstein condensate [7]. 
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Fig. 6. Contrast of Stiickelberg interference. Total transfer fraction ^ versus for = 1.000(l)7r. Values 
and error bars denote the mean and standard deviation of three consecutive measurements. The lines show the 
theoretical prediction for the fully integrated transfer fraction using the simple incoherent analytical model (blue 
solid line) and the extended model including the Stiickelberg correction factor (purple dotted line). The red dashed 
line includes a gaussian distribution of forces with a width of 60 Hz resulting from the harmonic trapping potential. 
The result from the numerical time-evolution of the trapped 2D system for N = 256 atoms is shown as the black 
dash-dotted line. 



the cloud size based on the trapping frequencies of approximately 30 pm one calculates a variation of 
the effective force over the cloud on the order of 60 Hz. We include this effect by assuming a gaussian 
distribution of forces over the entire sample, leading to different Stiickelberg oscillation frequencies and 
thus to a reduced contrast of the oscillations. The result is shown in Fig. 6, where the force averaged 
curve (red dashed line) is in close agreement with the experimental data. This picture is confirmed by the 
results of the numerical simulation of the trapped lattice system (black dash-dotted line), which already 
shows a comparable reduction in contrast of the oscillations for a strictly 2D system. While additional 
effects, such as the variation of the lattice depths along the third spatial axis z, are smaller than the effect 
identified above, they probably cause the remaining loss of contrast leading to no visible oscillations in the 
experimental data. Therefore it can be safely assumed that the simple incoherent version of the analytical 
model is sufficient for comparison with the experiment. 



4 Conclusion 



We have investigated the intricate dynamics of interband transitions arising from the presence of two 
Dirac points. Subsequent transfer through two linear band crossings leads to the appearance of a distinct 
feature: double peaks in the probability of finding atoms in the higher band after one Bloch oscillation, 
both as a function of sublattice energy offset as well as of quasi- momentum q^. We find good agreement 
with a theory without free parameters, based on a universal Hamiltonian describing the vicinity of the 
Dirac points [19,20]. The absence of Stiickelberg interference was successfully explained by large differences 
in the acquired Stiickelberg phases due to the inhomogeneity of the force applied to the atoms - an effect 
of the harmonic trapping potential. This explanation is confirmed by numerical simulations of a trapped 
2D lattice system. 

The momentum-resolved measurement of two sequential Landau- Zener transitions has been shown to 
be a sensitive probe of the band structure of our system. In this paper, the effect of varying the geometry 
of the lattice and the energy offset between sublattices was characterized in this way. In future work the 
infiuence of factors which are more challenging to model theoretically, such as impurities, interactions, 
finite size or different types of edges could be studied in the same manner. Furthermore, cold atoms 
in optical lattices may become a useful tool for investigating open questions in relativistic quantum 
mechanics using the linear dispersion relation surrounding a Dirac point. 
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A Numerical simulation 

We describe the optical lattice with the following tight-binding Hamiltonian [20], 



i^lattice = -h 4(^j Yl Yl 

{ij)er {ij)^b {ij)^9 



zeA 



ieB 



(5) 



where the hopping amplitudes of the red, blue and green bonds (r,b,g) are ^1,^2 and ts respectively, 
see Fig. _7. The lattice consists of two sites per unit cell. Aab is the staggered onsite energy on the A 
and B sublattices. The tight-binding parameters ti, ^2, ts and Aab are extracted from fits of the tight- 
binding Hamiltonian to the band structure calculated for the potential given in Eq. 1. The real-space 
Hamiltonian (Eq. 5) is Fourier transformed to momentum space to get a two-by-two matrix H\^. It is 
diagonalized by tjl^H\Jj\^ = E\^. When Aab = and ^2+^3 < 2ti, the tight-binding model features Dirac 
points [23]. 



Fig. 7. The tight binding model. Each unit cell contains two sites [A and B). The hopping amplitudes for 
the red, blue and green bonds are ti,t2 and ts respectively. For ts = the model reduces to a brick- wall lattice, 
which is topologically equivalent to the honeycomb lattice. 



In the experiment there is an additional harmonic trapping potential 



trap 



(6) 



where Xi^i/i are the spatial coordinates of the i-th site. They are measured in units of A/2. 7^^ 



iy) 



x{y) 



(A/2)^ are the strengths of the harmonic confinement along the x (y) direction. They are deter- 
mined from the intensities of the laser beams (see Appendix B). 

To simulate the Bloch-Zener oscillation experiment, we first solve the ground state \^) of N = 256 
spinless fermions of i^iattice + ^trap , then apply a linear gradient field infield = FX/2 yiUi to the system 
and evolve the wave function 



(7) 
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The time evolution is performed with the Lanczos algorithm [24] with 200 Lanczos vectors. The number 
of lattice sites (2 x 200^) is chosen such that the cloud does not touch the boundary during the Bloch os- 
cillation. The time step is O.OST^. At each time step, we measure the density matrix pij = {^{t)\cjCi\^{t)) 
of the system. 

To extract momentum distributions, we reshape pij into pjj, where /, J are indices for the unit cell, 
and a, b are sublattice indices. A Fourier transform with respect to /, J gives = pk- Applying the 
unitary transformation Ul-Pi^Ui^, the two diagonal elements become the band filling nk, lower and nk, upper- 



t=OTB t = 0.25TB t=0.5TB t = 0.75 Tb t = 1 Tb 




Fig. 8. Time-resolved Bloch oscillations. Comparison of the quasi-momentum distribution for Bloch oscil- 
lations resulting from a force pointing along x in the experiment (top) and in the numerical simulation of a 2D 
trapped system for N — 256 atoms (bottom). 

As a demonstration of the method, a comparison of one time-resolved Bloch oscillation cycle for a 
force pointing along x is shown in Fig. 8. As can be seen in the figure, the simulation captures the essential 
features of the experimental data. 



Numerical simulation parameters 

Fig. 3b. Left simulation: ^1,2,3/^ = [589, 969, 184] Hz, -fx.y/h = [0.09, 0.28] Hz, AAs/h = 20 Hz, FX/2h = 
63 Hz, right: ^1,2,3//^ = [620,525,131] Hz, -/^^y/h = [0.09,0.28] Hz, AAs/h = 17Hz, FX/2h = 63 Hz. 

Fig. 6, black dash-dotted curve. ^1^2,3 as obtained from a fit to the calculated band structure for the 
given lattice parameters, j^^y as deduced from the lattice parameters, see Appendix B. Aab/^ ~ 
20Hz, FA/2/i = 89Hz. ' ~ 

Fig. _8. V^/Er = 4.3(2), = 1.013(l)7r, corresponding to ti 2 3 = [602, 789, 166] Hz, 7^ y/h = [0.15, 0.59] Hz, 
Aab = 266 Hz. FX/2h = 89 Hz. 



B Lattice potential 

The phase (p (see Eq. 1) between the X and Y beams at the position of the atoms is stabilized interfero- 
metrically using a pair of additional beams detuned from each other and from X, X and Y. This results 
in a weak additional lattice along each axis of about OAEr. 

The laser beams used to create the lattice potential cause the atoms to experience a weak har- 
monic trap in all three spatial directions. For a lattice with V-^ x y/-^R ~ [4.0(2), 0.28(1), 1.8(1)], the 
trapping frequencies are given by uj^^y^zl'^^ — [17.6(1), 31.8(5), 32.7(5)] Hz, as calibrated from dipole 
oscillations of the cloud. In general, the trapping frequencies approximately scale as uj^ oc ^/Vy ^ Uy ex 
Y^V^+ {VxVy /Vx) and ujz oc y^^^^^^T^Oc^VT^^I^^^^^^^IV with respect to the intensities of the lattice 
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beams. The harmonic confinement causes the time it takes for the atomic cloud to reach the center of 
the B.Z. after a Bloch oscillation to be slightly larger than Tg. 



C Momentum-space integration of transfer fractions 

In the following we describe how the total transfer ^ for the entire atomic cloud in 3D is calculated 
from the momentum-dependent transfer £,d{Qx) following the procedure presented in [20]. We assume a 
semi-classical expression for the energy of the particles 



/ \ , % , Qz , ^ C 2 2, 2 2, 2 2\ ro\ 

'*^"2^^2^^2^^^ [m,oj,x + myujyy + m,uj,z ) , (8) 

X y 



where the effective masses rux^y^z are obtained by expanding the dispersion relation around quasi- 
momentum q = in the tight-binding regime and c^Ox^y^z are the trapping frequencies in the three different 
directions. At zero temperature the expression for the integrated transfer fraction is then given by 

Here we used the fact that the transfer only depends on and denote the Fermi wave vector qr = 
y/2m^/h and the chemical potential ja. For the semi-classical expression of the energy of the particles, 
the total atom number is related via /i = fuj{6Ny^^ , where uJ is the geometric mean of the trapping 
frequencies. 
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